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ABSTRACT 


A  theoretical  analysis  of  the  plane  flow  produced 
by  the  relative  normal  motion  of  two  parallel  plates  is 
presented.  Creeping  motion  of  a  Newtonian  liquid  is  con¬ 
sidered  so  that  a  linear  system  results.  The  classical 
solution  is  presented  and  shown  to  lead  to  rather  artificial 
boundary  conditions.  These  boundary  conditions  are  modified 
by  the  superposition  of  an  auxiliary  solution  which  is  ob¬ 
tained  by  expanding  the  stream  function  in  a  double  Fourier 
series.  For  large  plate  width  to  separation  ratio  the 
auxiliary  solution  is  found  only  to  affect  the  flow  in  the 
region  of  the  plate  edges. 
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NOMENCLATURE 


U 

A 


u 

e 


i 

ij 


u,v 


X,Y 

U,V 

f 

I 


coefficient  of  viscosity 
coefficient  of  viscosity 
density 

stress  tensor  referred  to  rectangular  coordinate  axes  ox^ 
stress  deviator  tensor 
thermodynamic  pressure 


velocity  vector 
,  /du  .  du 

2  +  ^ 


rate  of  deformation  tensor 


N  J  i" 

rectangular  cartesian  coordinates 

velocity  components  in  the  x  and  y  directions 

respectively 

nondimensional  rectangular  cartesian  coordinates 
nondimensional  velocity  components 
stream  function 

nondimensional  stream  function 
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CHAPTER  I 


INTRODUCTION 

This  thesis  considers  isothermal,  creeping,  laminar 
flow  of  an  incompressible  Newtonian  fluid.  Although  many 
applications  of  fluid  mechanics,  such  as  hydrodynamic  lub¬ 
rication,  involve  heat  dissipation,  the  assumption  of  iso¬ 
thermal  flow  eliminates  the  complication  of  variable  viscosity 
due  to  temperature  gradients  and  does  not  impose  a  great  re¬ 
striction  on  the  use  of  the  analysis.  Creeping  flow  is 
considered  and  this  results  in  linear  differential  equ¬ 
ations.  However  one  of  the  objects  of  this  thesis  is  to 
indicate  a  scheme  for  the  inclusion  of  inertia  effects. 

With  the  governing  equations  linear,  the  superposition 
principle,  which  is  often  used  in  classical  elasticity  prob¬ 
lems,  can  be  applied. 

Since  liquids  show  negligible  compressibility  effects, 
even  under  extreme  pressure,  an  incompressible  fluid  with 
ideal  Newtonian  properties  is  the  mathematical  model  upon 
which  the  analysis  of  this  thesis  is  based.  Many  similar 
problems  have  been  solved  subject  to  these  restrictions 
and  the  results  found  to  be  in  good  agreement  with  experi¬ 


mental  observation. 
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( 8 ) 

The  classical  solution  given  by  Michellv  '  for  creeping 
laminar  flow  between  parallel  plates  with  relative  normal 
motion  implies  a  parabolic  distribution  of  the  radial  vel¬ 
ocity  component „  This  solution  satisfies  the  Navier  Stokes 
equation  and  the  velocity  boundary  condition  at  the  plates 
but  gives  a  rather  artificial  stress  distribution  at  the 
edges  * 

For  large  radius  to  separation  ratios  a  principle 
analogous  to  St.  Venant ' s  principle  in  classical  elasticity 
can  be  applied  and  these  boundary  effects  neglected.  A 
more  satisfactory  solution  can  be  obtained  by  superimposing 
a  solution  on  the  parabolic  profile  and  thereby  satisfying 
the  natural  boundary  condition  of  constant  normal  stress 
and  zero  shear  stress  at  the  edges. 

To  avoid  the  complexity  of  cylindrical  coordinates 
the  plane  problem  is  considered.  The  origin  of  the  coor¬ 
dinate  system  and  the  symmetrical  movement  of  the  plates 
were  chosen  to  facilitate  a  possible  later  study  of  inertia 
effects.  For  the  quasistatic  problem  the  solution  is  not 
affected  if  a  velocity  of  translation  is  superimposed  so 
that  one  plate  is  stationary.  With  both  plates  moving  the 
origin  is  taken  at  the  centre  in  order  to  keep  the  frame 
of  reference  stationary. 

The  plane  problem  is  solved  using  a  Fourier  series 
expansion  of  the  stream  function.  The  governing  equations 
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necessitate  differentiating  this  series  four  times.  In 

( 3  4) 

two  papers  on  the  stability  of  viscous  flow  Goldstein  * 
shows  how  Fourier  series  can  be  used,  the  justification 
for  differentiation  depending  on  the  conditions  satisfied 
by  the  function  at  the  end  points  of  the  interval.  Green 
shows  how  double  Fourier  series  can  be  applied  to  boundary 
value  problems  that  involve  fourth  order  partial  differential 
equations.  The  use  of  a  double  series  avoids  some  of  the 
algebraic  difficulties  of  a  single  series. 

In  order  to  satisfy  the  boundary  conditions  of  the 
problem  directly,  the  expansion  of  the  normal  velocity  com¬ 
ponent  at  the  plates  in  a  single  series  is  required  so  that 
it  may  be  incorporated  in  the  solution.  A  few  terms  of  the 
resulting  infinite  series  shows  the  Gibbs  Phenomenon  which 
is  due  to  non  uniform  convergence  at  x  =  +  a.  Even  after 
many  terms  are  used  the  velocity  profiles  do  not  show  con¬ 
servation  of  volume  and  an  overshoot  develops  at  one  point 
on  the  profile.  The  development  is  due  to  the  Gibbs  Phen¬ 
omenon,  which  arises  in  the  expansion  of  discontinuous 
functions „ 

Specified  stress  boundary  conditions  at  x  =  +  a 
can  be  satisfied  by  superimposing  an  auxiliary  solution 
on  the  parabolic  solution.  The  auxiliary  solution  must 
satisfy  the  boundary  conditions  of  zero  velocity  at  y  =  +  b. 
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CHAPTER  II 


GOVERNING  EQUATIONS  OF  VISCOUS  FLOW 


The  governing  equations  for  laminar  flow  of  an  in¬ 
compressible  Newtonian  viscous  fluid  are  the  continuity 

equation  _ 

V  •  U  =  0 

where  U  is  the  velocity  and  the  Navier  Stokes  equation 

dU  _  ,  -  ,  2- 

p  dt  _vp  +  PX  +  H  V  U 

where  X  is  the  body  force  per  unit  mass. 

These  equations  applied  to  the  two-dimensional 

j 

finite  region  bounded  by  two  parallel  plates  (Fig.  2.1) 
in  the  absence  of  body  forces  are 


du  ,  dv 

^  +  ^ 


0 


/du  ,  du  ,  du\ 

p  +  u  ^  +  V^  J 

/du  ,  dv  ,  dv\ 

p  [it  +  u  ^  +  v 


(2.1) 


(2.2) 


(2.2b) 


(2.3) 
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Fig.  2.1 

The  equations  (2.3)  can  be  reduced  to  non-dimension¬ 
al  form  with  the  use  of  a  suitable  characteristic  velocity 
and  length  for  the  problem.  If  the  plates  are  allowed 
to  separate  with  constant  relative  velocity  V,  their 
width  being  2a,  the  characteristic  velocity  and  length 
become  W  and  a. 

The  non-dimensional  variables  become 


t 


tw 

a  9 


y 


.  2L 


a 


u 


v 


Since  a  and  b  are  considered  in  this  problem  to 
be  of  the  same  order  all  of  the  non-dimensional  variables 
are  0(1). 


■ 


' 


-noxe.  ti:fo— non  od  b  j  i  nto  (£•£  snox:!6'rp9  or  r 


.it  j  srfd  xo.t  dipnoi  n- 


6 


Substituting  these  quantities  equations  (2.3)  become 


du 

35? 


+ 


U 


,  —  dv  — 
+  u  3=  +  V 
ox 


(2.5) 


where 


M- 


R  the  Reynolds  number. 


For  creeping  motion  the  Reynolds  number  approaches 
zero  and  since  the  terms  in  brackets  are  0(1)  the 
material  derivative  can  be  neglected.  This  is  a  form 
of  Newtonian  flow  in  which  the  pressure  and  viscous  force 
are  of  the  same  order  and  the  inertia  term  is  neglected. 

The  approximation  involved  by  neglecting  the  inertia 
term  can  be  justified  if  the  solution  obtained  by  making 
the  approximation  satisfies  the  approximation  used  in 
obtaining  it.  This  will  be  investigated  in  the  conclud¬ 
ing  chapter  and  the  simple  solution  obtained  shown  to 
satisfy  the  approximation. 
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For  creeping  motion  where  the  material  derivative 
is  neglected  equation  (2.2)  becomes 

-  vp  +  pX  +  |x  V2U  =  0  .  (2.8) 

Considering  only  body  forces  which  are  derivable  from 
a  scalar  potential  <t>  equation  (2.8)  becomes 

V(p  -  p<t>)  =  ]±  V2U  (2.9) 

and  for  plane  flow 

ar  (p  _  p*] 

.  ^  (P  -  pO) 

From  equation  (2„9) 

curl  V2U  =  0 

or 

_  d3u  _  d3u  +  d3y  +  d3v 
dx2dy  dy3  dx3  dxdy2 


=  H 


'd  2_u  +  d2u 


>x 


dy‘ 


=  H 


^d2v  +  aV 


ix 


8y^ 


(2.10) 


0 


(2.11) 
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The  continuity  equation  implies  the  existence  of  a 
stream  function  such  that 


Substituting  equations  (2.12)  in  (2.11)  gives 


2 


dx2dy2 


0 


> 


or 


0  . 


This  is  the  differential  equation  for  the  stream  funct 


for  creeping  flow  of  an  incompressible  Newtonian  fluid. 


(2.12) 


(2.13) 
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CHAPTER  III 


USE  OF  FOURIER  SERIES 


The  method  used  in  this  thesis  involves  the  differen¬ 
tiation  of  a  Fourier  series  four  times.  In  many  applications 
of  Fourier  series  the  justification  for  differentiation  is 

often  overlooked.  In  two  papers  on  the  stability  of  viscous 

( 3  4) 

flow  Goldstein  s  showed  how  Fourier  series  can  be  used, 
the  justification  for  differentiating  depending  on  the  con¬ 
ditions  satisfied  by  the  function  at  the  end  points  of  the 

( 5) 

interval.  Green  extended  Goldstein's  work  to  double 
Fourier  series. 

If  F(x,y)  is  defined  and  bounded  in  the  range  (-tt  —  x  ir, 
— 7T  ~  y  ^  tt)  and  if  F(x,y)  satisfies  the  conditions: 

(a)  F(x,y)  has  only  a  finite  number  of  maxima  and  minima 

(b)  F(x,y)  has  only  a  finite  number  of  finite  discontinuities, 
it  may  be  expanded  as  a  double  Fourier  series.  These  con¬ 
ditions  are  sufficient  and  are  valid  for  the  functions  which 
arise  in  this  application  to  fluid  flow. 

If  F(x,y)  satisfies  the  Dirichlet  conditions  and  if, 
outside  the  range  (-tt  ^  x  ^  t r,  -tt  ^  y  ^  it)  F(x,y)  is  periodic 
in  the  sense  that 


F(x  +  2tt,  y)  =  F  (x,y) 
F  (x,  y  +  2tt)  =  F  (x,y) 


.bsdool'xsvo  naiio 


t  *  so  ii  i  h  n  b  i,  uJ  >.n  ;  ;  '  (y  -  S  C 


■nirtim  bfTjs  Bmixwfl  io  i9o[mun  s  ;  rnil  6  yin  :>  a /rf  (y  )  n 

' ],tr. 

.391193  isi  ’UO^  o  djjob  £  as  bebn.  ••  3  5C  r:  3 


•  no,  o  [ui  srfd  loi  biLsv  ois  bne  dnsioiirfcifa  ai£  ?  icx^ib 


(\  «:.)  [  *  (Y  £  +  :)' 


10 


then  F  can  be  represented  by  a  double  trigonometric  series. 

For  example  a  function  F(x,y)  that  is  an  odd  function  of  x 
and  y  and  has  half  range  symmetry  in  x  and  y  can  be  represented 
by  oo  oo 


F 


(*,y) 


A  sin  rx  sin  sy 
rs  2 


r  =  1, 3  .  .  .  s  =  1,3... 


where  the  A  are  defined  by 
rs  J 


rs 


16 

2 
t r 


7T  P  7T 

2  /  2  F(x,y)  sin  rx  sin  sy  dxdy 


o 


For  a  fixed  y  the  series  converges  to  the  sum 


2  ,F(x  +  0,  y)  +  F(x  -  0,  y) 

L.  _ 

where  F(x  +  0,  y)  is  the  value  of  F  as  x  is  approached  from 
the  larger  value  and  F(x  -  0,  y)  is  the  value  of  F  as  x 
is  approached  from  the  smaller  value.  Similarily  for  a 
fixed  x  the  series  converges  to 

F(x,  y  +  0)  +  F  (x,  y  -  0) 

If  F(x,y)  is  continuous  at  the  point  this  reduces  to  F(x,y). 
Also  due  to  the  symmetry  of  the  problem  only  odd  values  of 
r  and  s  will  appear  and  will  not  be  shown  explicitly. 

If  F(x,y)  is  a  function  which  can  be  expanded  in  a 
double  Fourier  sine-sine  series  in  the  region  —  x  ^  , 

-  ^  rE-  y  <  ,  and  if  its  partial  derivative  dF/dx  can  be 

expanded  in  a  cosine-sine  series,  the  coefficients  in  each 


- 
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case  being  formed  by  the  usual  rule,  then  the  coefficients 
in  the  second  series  can  be  related  to  those  in  the  first 
series  by  a  partial  integration.  Similar  remarks  are  valid 
for  sine-cosine,  cosine-cosine,  and  cosine-sine  series. 

Thus  if 

F  (x,y)  =  Y_Y_  Ars  sin  rx  sin  sy  (‘  f  ~  x  ~  f »  ~  f  —  y  ^  f ) 


then 

5F{x,y) 

OX 


cos  rx  sin 


sy 


Similarly  if 
F (x,y) 


cos  rx  sin 


sy 


—  <•  x  <1  - 
2  *  2  * 


and 

dF (x.y) 

OX 


sin  rx  sin 


7T 

2> 


The  B  are  found  by  the  calculation 
rs 


B 


rs 


_  16 


t r 


jr  t r 

2 


2  sin  rx  sin  sy  dxdy 


7 r 


16 

2 

7T 


F  sin  rx 


7T  7T 

2  r  2 


J 


Fr  cos  rxdx  )  sin  sydy 


a  j.  16 

=  -r  A  +  —r 

rs  2 
7 r 


7T  r-1 

2  F^~,  y J  (-1)  2  sin  sydy 


. 
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Thus 


dF(x.y) 

OX 


r-1 

r  A  +  (-!)  2  d 

rs  i 


sin  rs  sm  sy 


2T  ^  x  — 

9  —  —  9  » 


where  d 


16 

2 
t r 


7T 

2 


F  f-  y. 


sin  sydy 


The  coefficients  in  the  double  Fourier  series  expansions 
of  F(x,y)  and  dF(x,y)/dy  can  be  found  similarly. 


CHAPTER  IV 


METHOD  OF  SOLUTION 


4.1  Preliminary 


The  classical  solution  given  by  Michell  for  creeping 
laminar  flow  between  parallel  plates  with  relative  normal 
motion  is  not  completely  satisfactory.  Although  the  vel¬ 
ocity  field  satisfies  the  Navier  Stokes  equations  it  re¬ 
sults  in  a  rather  artificial  stress  condition  at  the  edges. 
Incompressible  flow  is  considered  since  this  is  a  realistic 
assumption  for  a  liquid.  Inertia  forces  are  neglected  but 
body  forces,  obtained  as  the  gradient  of  a  scalar  potential 
function  can  be  incorporated. 

The  problem  considered  is  the  modification  of  the  plane 
flow  analogue  of  Michell 's  solution  by  the  superposition 
of  an  auxiliary  problem  in  order  to  satisfy  the  boundary 
condition  of  constant  normal  stress  and  zero  shear  stress 
at  the  edges. 

y 


FIG.  4.1  COORDINATE  SYSTEM 
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The  origin  0  of  the  coordinate  system  is  taken  as 
shown  in  figure  4.1  to  maintain  a  stationary  coordinate 
system.  The  plates  move  with  symmetrical  normal  motion, 
relative  to  the  coordinate  system,  either  apart  or  together. 
For  the  purposes  of  the  analysis  they  are  considered  to  move 
apart  but  the  reverse  effect  involves  only  a  trivial  change 
of  sign.  Although  the  quasistatic  case  is  not  affected, 
except  for  a  superimposed  velocity  of  translation,  by  mov¬ 
ing  only  one  plate,  the  symmetrical  motion  is  desirable  for 
possible  future  consideration  of  inertia  effects  with  the 
separation  velocity  a  function  of  time. 

The  governing  equations  are 


du  ,  dv 

^  +  57 


0 


d 

3^ 


(p 


d 

dy 


(2.2b) 


(2.10) 


and  the  boundary  conditions  are 


u  =  0 
W 


at 

at 

at 


y  =  +  b 
y  =  b 
y  =  -b. 


(4.1  a ,b , c) 


The  imposed  boundary  condition  of  no  slip  at  a  solid 
boundary  has  been  a  subject  of  much  discussion.  Under  most 
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circumstances  this  condition  has  been  experimentally  veri- 
fied  although  it  cannot  be  proved  from  hydrodynamical  con¬ 
siderations  ,  This  boundary  condition  may  not  be  realistic 
for  extreme  conditions  of  rarefied  gases  and  very  high  vel¬ 
ocity  where  molecular  effects,  which  are  not  included  in 
continuum  theory,  are  significant.  In  hydrodynamical  lub¬ 
rication  the  no  slip  boundary  condition  is  assumed  to  hold 
except  where  the  film  thickness  becomes  comparable  with  the 
mean  free  path  of  the  molecules.  This  microscopic  effect 
is  also  excluded  from  consideration, 

4,2  Elementary  Solution 

The  no  slip  boundary  conditions  (4,1a)  are  satisfied 
by  the  expression 

u  =  cx(y2  -  b2)  (4,2) 

Substitution  of  (4,2)  in  the  continuity  equation  (2,2a)  gives 

!  2  ,  2v  ,  dv  _ 

c(y  -  b  )  +  ^  -  0  , 

Consequently 

v3  2 

v  =  -c  +  cb  y  +  f  (x) 

and  f(x)  =  0  since  v  is  independent  of  x  at  y  =  +  b  and 

v  =  0  at  y  =  0 . 


w,  -i  yXXfc^n  nx^3qxe  n9^cf  serf  aoi^x6n  ?o  sin -j  eeon  s  t^nu/oiio 

-X =>v  rfpirf  v  lev  brie  89SBp  bsx^dSBtx  3to  s  oxiibnoo  jffisiixe  to :* 

'  ■ 


no±:tuXo8  YiB^nsmaXa 


■ 


anonixbnoo  yisbnuocf  qxls  on  3ffT 
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From  (4.1b) 


c 

3W 

4b3 

therefore 

3W 

u  =  —  3 
4bJ 

x(y2  - 

b2) 

and 

3W 

v  =  — 

4d 

y(b2  - 

2 

(4,3) 


(4,4) 


Equations  (4.3)  and  (4.4)  satisfy  the  Navier  Stokes 

2- 


equations  since  curl  (v“U)  =  0.  Substitution  of  (4.3)  and 
(4.4)  in  the  Navier  Stokes  equations  (2.10)  gives 

(ip  _  3  Wnx 

I  II  iTTi  av  f  ■■  i  n  ■ 


dx 

a 


2  b3 


n  - 


2 

2 


M 

b3 


Substitution  of  equations  (4,5)  and  (4.6)  in 

dp  =  dx  +  dy 

gives  the  exact  differential 


dp 


and  integration  gives 


3Wg 

2b3 


(x  dx  -  y  dy) 


3Wll  /  2  2  ,  _  \ 

p  -  (x  -  y  +  c  ) 

4b 


At  x  =  +  a  ,  y  =  +  b  take  p  =  0 


(4.5) 


(4.6) 


2  ,  2 
a  +  b 


then 


<:>•*  JG-  , 


17 


and 


-  ^  (X2  -  y2  -  a2  +  b2: 

4b 


A 


and  at  x  =  +  a 


3Wjljl  .  2  ,2, 

— o  (y  -  b  , 

4b 


The  normal  stress  component  in  the  x-direction  can  be 

j 

written  as 

a  =  s  -  p 

X  X  c 


and  from  equation  (2„5] 


x 


0  du 
2ti  ^ 


The  partial  derivative  with  respect  to  y  of  the  normal 
stress  is 


do 


x 


X 


408: 


and 


ds 

x 

civ 


2(1 


A2 
o  u 


'4  0  9' 


Substituting  (4,3)  in  (4,9)  gives 


x  


dy 


3uWy 

b3 


and  substituting  this  equation  and  (4.6)  in  (4,8)  gives 


do 


x  


"Sy 


9u.Wy 

2b3 


at  x  =  +  a  . 


(4.10) 


bnfi 


i  !  ,  ;\:v 


rfitw  sviJEVXisb  s/fT 
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The  shear  stress  is 


xy 


/du  ,  dv' 

*  1 55  +  ^ 


and  using  (4.3)  and  (4.4)  this  gives 


xy 


3  llWxv 

2b3 


and  at  x  =  a 


xy 


2b' 


Way 


(4.11) 


(4.12) 


The  distribution  of  stress  indicated  by  the  elementary 
solution  is  somewhat  artificial.  A  more  realistic  boundary 
condition  at  x  =  +  a  would  be  constant  normal  and  zero 
shear  stress.  Since  the  governing  equations  are  linear 
an  auxiliary  solution  can  be  superimposed  on  this  solution 
and  thereby  satisfy  the  conditions  of  constant  normal 
stress  and  zero  shear  stress  at  the  edges. 


4.3  Auxiliary  Solution 


The  range  (-a  x  ^ 
to  the  range  ( -  x  ' 

_  2aX 
x  =  — 

7T 

Using  the  dimensionless 


a,  -b  v  ^  b)  can  be  converted 
~  ~  ^  Y  ^  )  with  the  substitutions 

,  y  =  •—  (4.13) 

5>  ^  7r 

stream  function 


¥ 


Jl 

aW 


4.14) 


the  biharmonic  equation  (2.13)  becomes 


■ 


■ 


' 


' 
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a4y 

ax4 


2a‘ 


a4f 


ax2dx2 


+ 


a 


a4f 

ay4 


0 


(4.15) 


where  a  =  ^  is  the  ratio  of  the  plate  width  to  separation 
distance. 

Assume  that  Y  can  be  expanded  in  a  double  Fourier  sine 
series  as  follows 


¥ 


II 


A  sin  rX  sin  sY 
rs 


JL  <z  y  JL 
2  —  ~  2 

(4.16) 

JL  <  v  — 

2-2 


where  r  and  s  are  odd  positive  integers  because  of  the  sym¬ 
metry  of  the  problem 


3y  ' 


rs 


sin  rX  cos  sY 


d2y 

dY2 


2 

s 


A 

rs 


sin  rX  sin  sY 


3 

s 


A 

rs 


sin  rX  cos  sY 


s-1 


sin  rX 


sin  sY 


JL  ^  ^  JT 

2  ^  2 


JL  v-  — * JL 
2  2 


JL  ^  y  ^  JL 
2  ~~  2 


JL  y  <=-  — 
2  ~  2 


JL  v  <=z  JL 
2  —  A  —  2 


JL  y  <  JL 
2  2 


-  %■*=  x-^  % 

(4.17a) 

JL  <  y  JL 
2  X  2 


awollo^  36  esn  >8 


. 
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by  ... 

3x  " 


r  A  cos  rX  sin  sY 
rs 


-  ^  ^  X  <C  4 


7T 

2 


2  -  *  —  - 


7T 

2 


a2f 

ax2 


r-1 

■r2  A  +  (-1)  2  c  j  sin  rX 
rs  s  / 


7T 

2 


x-f 


sin  sY 


jr 

2 


Y  <  — 

~  2 


d3y 

ax3' 


r-1 

-r 3  A^g  +  r(-l)  2  cg J  cos  rX 


7T 

2 


7T 

2 


sin  sY 


_  JL  y  <c.  jl 
2  — - 


7T 

2 


rs 


r-1 

r2(-l)  2 


r-1 


cs  +  (_1) 


7T 

2 


4  ^  X  1= 


z 

2 


Z 

2 

Z 

2 


a4y 

dx2dY2 


r-1 

22_  2  ,  2 
r  s  A  -sc  (-1) 
rs  s 


Z 

2 


—  X 


z 

2 


sin  rX  sin  sY 


-  2  ~  Y  — 


Z 

2 


where 


b  = 


c  = 


d  = 


i6  r 

z 

d3¥ 

(X 

TT2  I 

2 

dY3 

o 

i6  r 

Z 

(Z 

v2> 

2  / 

2 

^X 

T  o 

m  r 

Z 

2 

d3¥ 

(Z 

v2 

^  jo 

dx3 

t r 


)  sin  sY  dY 


dY 


(4.17b) 

(4.17c) 

(4.17d) 

( 4 , 1 7  e ) 


( 4 . 17f ) 


. 
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Substituting  equations  (4.17a,b,c)  in 


+  Ct2s2)  A  +  a4  (-1)  2  b  -  (-1)  2 

XT  b  Jl 


c  -  d 
s  s 


=  0 


(4.15)  gives 
(r2  +  2a2s2) 


s-1 


r-1 


-a4  (-1)  2  \  +  (-D  2 


,2,^22 

(r  +  2a  s  c  -  d 

s  s-J 


rs 


,  2  ,  2  2,  2 

(r  +  a  s  ) 


(4.18) 


This  expression  is  valid  for  the  plane  flow  of  a 
liquid. 

The  boundary  conditions  for  the  auxiliary  problem  ares 


U  =  -  =  0  at  Y  =  +  ^ 

w  -  2 

V  =  -  =  0  at  Y  = 

w  -  2 


(4 . 18a ,b) 


and  making  the  stress  conditions  for  the  auxiliary  solution 
the  negative  of  the  stress  conditions  for  the  elementary 
solution  equations  (4.10)  and  (4.12)  yield 


-18ll  WY 


at  X  =  +  f 


-3aLL  WY 

Trb 


at  X-  f 


(4 . 18c ,  d) 


The  boundary  condition  (4.18a)  is  satisfied  by  the 
elementary  solution  and  use  of  (4.18b)  and 

V  =  -  ■—  r  Ars  cos  rX  s-*-n  sY  -  ^  <  X 


-  - 


' 


> 

22 


gives 


r 


S”1 


s 


cos  rX  “  0 


which  is  true  for  all  r 


Therefore 


s-1 


A  ( 
rs  ' 


0 


Equation  (4.17e)  yields 

2  r— i 

V  =  “  “8"  /  Cs  sin  sY  at  X  =  +  ” 

s 

and  the  use  of  (4018b)  gives 


s-1 

^  cg  (-1)  2  =  0  at  X  =  +  | 

s 

The  boundary  condition  (4018b)  is  satisfied  by  (4e19)  and 


(4 0 19a)  0 

Substituting  the  components  of  velocity  derived  from 
the  stream  function  (2.12)  into  equation  (408)  gives 


8ax 

~5y 


7 2  U.W  f 1  d  3f 

4b  Va2  *X3 


b3'¥ 

dxdY2 


Equations  (4.17e,f)  yield 


3x 


w 

4 


c  sin  sY 
s 


at  X 


7 r 


d3f 

Sx3 


jr 

4 


d  sin  sY 
s 


at  X  =  ~ 


TT 

2 


(4,19) 


(4„ 19a) 


(4 „ 20) 
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Hence  C 4  0  20)  becomes  at  the  edges 


da 

X  _ 

dY 


3 

Teb 


'■a 


s  ,  2  \  __ 

2  +  s  cs  j  sm  sY 


(4.2i; 


Expanding  the  boundary  condition  (4,18)  in  a  Fourier  sine 


series  gives 

WY 

2, 

7 r  b 


•72u  w 

3, 

7 r  b 


s~l 

2 


sin  sY 


(4,22) 


Equating  (4,21)  and  (4,22)  yields 


9  2 

d  +  as  c 
s  s 


1152  Q  Jtli 


s-1 

2 


7 r 


6 


=  0 


Eliminating  d  from  equation  (4,18) 

s 


rs 


s-1 

4 ,  .  >  2  , 

-a  (-I)  b 

j 

,2,2  2>2 

(r  +  a  s  ) 


r-1 

2 


s-1 


+ 


o  90  9  1152  (-1)  ‘ 

(r  +  3a  s  )c  -  a  6  2 

s  7 r  s 


,  2,2  2,2 

(r  +  a  s  ) 


(4,23) 


Similarly  the  shear  stress  condition  (4,11)  becomes 


tt2ll  W  /l  d2¥ 


xy 


4b 


dx2 


a 


dV 

dY2 


Substitution  gives 


=  7T  LL  W 
xy  4b 


-r2A 


rs 


a 


r-1 

2 


a 


c  -  a2s2A 
s  rs 


r-1 


(-1) 


sin  sY 


(4.24) 


' 


■ 
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Expanding  the  boundary  condition  (4018d)  gives 


■3aLi  WY 
Trb 


■12au  W 
2 

t r  b 


■1) 


s--l 

2 


sin  sY 


(4.25) 


Equating  (4C24)  and  (4.25)  yields 


r-1 


>-«  2  <r2  -  «V,  »r.  -  H  «2 

7 r  s 


s-1 

2 


=  0 


(4 0  26) 


Substitution  of  (4023)  gives 


r+s-2  _4 ,  2  _2  2V 

—  v  -  a  (r  -as  )b 

(_D  z  — - -  — - - — -£  + 

'  j  ,  2  .  2  2.  2 

(r  +  a  s  ) 


,24 
4a  s  c 


,  2  ,  2  2,  2 

(r  +  a  s  ) 


48* 

4 
7 r 


s-1 


idtl 


24  V1  (r2  -  a2s2)  ' 

2  /  ,  /  2,2  2,2 

7 r  ^  (r  +  a  s  )  ' 


(4.27) 


Substitution  of  (4.23)  in  (4.18)  gives 

(-1) 


a4b 


r+s  -  2 
2 


2  _i_  2 , 

r  +  3a  s  )  c 


,  2  ,22x2 

(r  +  a  s  ) 


s  


,  2  ,  2  2*  2 
(r  +  a  s  ) 


(-1) 


r-1 

2 


1156  2 

“6”  a 
7 r 


2  /  2  2  2,2 
s  (r  -as) 


(4.28) 


and 


s-1 

cs  (-1)  2  =  0 


(4.19a) 


Equations  (4.27) ,  (4.28)  and  (4.19a)  are  sufficient  to  de¬ 


termine  the  sets  of  b  and  c  from  which  the  set  of  coefficients 

r  s 


' 
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A  can  be  determined, 
rs 

The  effect  of  the  auxiliary  solution  on  the  total 
force  required  to  maintain  the  constant  relative  normal 
motion  of  the  two  plates  can  be  found  as  follows.  Since 
the  plates  represent  a  plane  of  maximum  shear  the  normal 
stresses  are  equal  and  the  stress  deviator  is  zero. 
Therefore 

a  =  a  =  ~p 
x  y 

The  Navier  Stokes  equations  (2,10)  become  on  substi- 
tut ion 


2 


+  a 


(4  e 29) 


Substitution  of  the  appropriate  derivatives  of  (4,17)  and 
setting  Y  =  ~  gives 


0 


at  *  "  | 


Since  U  =  V  =  0 


at  Y  =  " 


du 

c)X 


0  at  Y  =  - 


and  equations  (4,29)  yield 


.ous  a  ■  vob  srf:*  bn b  X6d  > s  a i  ceeasn  ?a 

' 


»  *  a 


Substituting  in  (4.29)  and  integrating  yields 


P 


-7 r3  a3u  W  V1 

16b  L 

r 


b 

cos  rX 
r 


at  Y  = 


which  satisfies  the  condition  of  zero  pressure  at 


roM 


CHAPTER  V 


CONCLUSION 


5„1  Summary  of  Results 

The  equations  (4„19),  (4e27)  and  ( 4  0  28)  were  set  up 

in  a  program  to  be  run  on  the  IBM  7040  computer .  This  in¬ 
volved  finding  the  coefficients  of  the  unknowns  b  and  c 

10  s 

and  arranging  them  in  a  suitable  form  for  solution  of  the 
simultaneous  equations „  The  program  was  written  in  double 
precision  to  obtain  sufficient  accuracy  when  a  large  num¬ 
ber  of  equations  were  used0  The  subroutine  Jordan,  which 
is  available  in  the  Department  of  Computing  Science  was 
used  to  solve  the  equations,. 

The  problem  was  first  solved  satisfying  only  the  nor¬ 
mal  stress  and  normal  velocity  conditions  at  the  edges 
(equations  (4019)  and  (4„28))0  A  possible  condition  which 
satisfied  the  normal  stress  condition  at  the  edges  was  to 
assume  all  the  cg  coefficients  to  be  zero  and  solve  equation 

(4„28)  for  the  coefficients  b  „  The  shear  stress  condition 

s 

(4027)  was  then  not  satisfied  and  the  shear  stress  for 
the  combined  problem  remained  approximately  equal  to  the 
condition  for  the  elementary  solution „ 

The  effect  of  adding  the  auxiliary  solution  is  to  in¬ 
crease  the  X  velocity  component  in  the  neighborhood  of 
the  midpoint  of  the  two  plates  and  decrease  the  velocity 
near  the  plates „  This  effect  dies  out  rapidly  for  increased 
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FIGURE  5.1  X  VELOCITY  PROFILES  FOR  PLATE  WIDTH/SEPARATION 
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FIGURE  5 „ 2  Y  VELOCITY  PROFILES  FOR  WIDTH/SEPARATION 
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FIGURE  5.4  Y  VELOCITY  FOR  WIDTH/SEPARATION  *  3  WITH  COEFFICIENTS 
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FIGURE  5 „ 6  Y  VELOCITY  FOR  WIDTH  SEPARATION  =  5  WITH  COEFFICIENTS 
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ratio  of  plate  width  to  separation.  For  larger  plate  width 
to  separation  the  addition  of  the  auxiliary  solution  pro¬ 
duced  only  small  changes  in  the  classical  solution  and 
these  changes  were  negligible  at  distances  from  the  edge 
greater  than  the  separation.  The  maximum  effect  was  approx¬ 
imately  a  ten  percent  change  in  the  velocity  distribution. 

In  order  to  satisfy  the  shear  stress  condition  (4.27) 
equation  (4,19)  was  substituted  for  the  last  of  equations 
(4.27)  and  with  equation  (4.28)  the  resulting  matrix  was 
solved  for  the  c  and  b  .  For  width  to  separation  equal 
to  one  the  velocity  profiles  in  the  interior  region  agreed 

with  the  solution  with  the  c  zero.  At  the  edges  the  vel- 

s 

ocities  did  not  form  a  smooth  curve  with  the  c  not  equal 

s 

to  zero  and  is  possibly  caused  by  the  slow  convergence 
of  the  solution.  The  results  were  also  unsatisfactory  for 
larger  width  to  separation  ratios. 


5.2  Discussion 


The  problem  solved  in  this  thesis  has  completely 
ignored  inertia  effects.  It  is  now  shown  that  the  vel¬ 
ocity  components  found  by  making  this  approximation  satisfy 
this  approximation. 

The  x  component  of  velocity  for  the  elementary  solu¬ 
tion  is 
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x (y  -  b  ) 
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Substituting  this  component  in  the  material  derivative 
and  the  viscous  term  and  letting 
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it  is  found  that 


du 


UV  u 


_  £Wb  /  3  (y  2-l) 

n  V  4 


1  +  3  (y 

2  8 


(y2-i) 


A  similar  result  can  be  found  for  the  y  velocity  and 
it  can  be  seen,  since  the  quantities  in  brackets  are  0(1), 
that  the  neglect  of  the  material  derivative  is  justified 
for  Reynolds  number  approaching  zero. 

In  the  one  dimensional  problem  for  radial  flow  between 

( 7) 

stationary  plates  J.L.  Livesay  has  shown  that  the  inertia 
effects  can  be  significant  at  relatively  low  velocity.  Un¬ 
like  the  problem  considered  the  characteristic  length  was 
taken  as  the  separation  between  the  plates  which  was  small 
compared  with  the  radius  of  the  plates.  A  method  such  as 
that  suggested  by  Karman  in  this  momentum-equation  approach 
to  boundary  layer  theory  can  be  used  to  incorporate  inertia 
effects.  A  velocity  profile  such  as  the  one  developed  in 
Chapter  IV  can  be  assumed  and  the  terms  of  the  equation  of 
motion  integrated  across  the  separation  between  the  plates. 
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APPENDIX 


COMPUTER  PROGRAMS  AND  OUTPUT  DATA 

FOR  VELOCITY  FIELDS 

The  calculations  for  the  velocity  fields  correspond¬ 
ing  to  various  plate  width  to  separation  ratios  was  greatly 
facilitated  by  the  use  of  an  IBM  7040  digital  computer. 

The  solution  of  the  large  number  of  equations  required  was 
well  suited  as  a  computer  problem.  The  Fortran  source 
programs  were  designed  for  a  general  number  of  equations 
so  that  the  effect  of  varying  the  number  of  equations 
could  be  examined. 

The  coefficients  of  the  unknowns  b  and  c  in  equations 

10  s 

(4.27),  (4.28)  and  (4.19)  were  arranged  in  matrix  from  suit¬ 
able  for  the  solution  of  simultaneous  equations.  The  num¬ 
ber  of  equations  solved  was  taken  equal  to  2L  allowing 
the  determination  of  L  of  the  unknowns  b^  and  L  of  the 

unknowns  c  .  These  equations  were  composed  of  L  equations 
s 

satisfying  the  normal  velocity  conditions  (4.28) ,  L-l 

equations  satisfying  the  shear  stress  equations  and  one 

equation  satisfying  the  normal  velocity  condition  at  the 

edges.  The  subroutine  Jordan  available  in  the  Department 

of  Computing  Science  library  was  used  to  solve  for  the 

unknowns.  The  coefficients  A  were  then  determined  from 

rs 

equation  (4.23)  and  the  accuracy  of  the  derivation  and  the 


' 

■ 
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precision  of  the  solution  were  then  checked  by  substituting 
the  A  into  equations  (4.18)  and  (4.23) .  The  program  and 

3T  S 

subroutine  were  written  in  double  precision  and  only  when 
a  large  number  of  equations  were  used  was  it  found  that 
the  number  of  significant  figures  declined  below  a  number 
acceptable  for  the  problem. 
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FORTRAN  SOURCE  FTATFMFNT  FOR  T  HF  VELOCITY  FIELD 


DOUBLE  PRFCTSTON  A(50*51>*  P T * P , R , S * SUMR * SUMRR , SUMS  9 SUMSS * P T ♦ PTP « 
IPR(lO)  *  P  R  P  (  1 0  )  *X*Y*W*Z*UR( 10) ,US(10) 10) *VR( 10) *VS(1 O) ,V( 10) * 

1  UP  no  )  *  VP  (  1  0  )  *  C  (  2  5  )  ,  R I  ( )  »RT I ( ?B  )  ♦  A  T  (OF)  ,  ARC ( CO ,F0 )  ,CI  (  ?5  )  *R ( ?5  ) 
INTFGFR  NJ*K*L*M*N*RHO*F,F*H 
PI=?.141 59 265 2 5 807^2 
DO  1 00  RHO= 1  >5,2 
8  P=PHO 
L  =  1 0 

M  =  7*(_  +  l 

DO?  T  =  1  *  L 
Y  =  L+I 
DO"*  J=1  *L 
F= I +J-? 

R  =  2*  J— 1 
S  =  2  *  I  —  1 

A ( T ♦J)=(-1D0)**F*(P**?-(P*S)**2) /( (R**?+(D*S)**?)**?)*P**4 

A(L*J)=ODO 

R=?*T-1 

S=2* J-l 

N  =  L+J 

A  (  K  ,N  )  =-  (-1  00  )  **f*  (  p**?  +  'anO*  (d*C)**7)  /  (  (p**?+(p*c)**?)r*?) 

A ( T ,N )=0DO 
F  =  J-l 

A ( L  *N )  =  ( -IDO ) **F 
7  A(K,J)=ODO 
!  IMP  =  ODO 
C(  |MRR  =  OD0 

C|  |mc  =  ODO 

Cl  JMCCaODO 
DO  15  J  =  1  ,L 
S=?*I-1 
R  =  2  *  J  —  1 

SUMRR=SUMRR+(R**2- (P*S )  **? )  /  (  (R**?+  (P*s  )  **?  )  **7  ) 

SUMR=SUMR+( (R**?-(P*S )**?)*( R**2+2 DO* (P*S)**2 ) ) / ( ( R**?  + ( P*S ) **? ) ** 
1  ?  ) 

P  =  2*  T-l 
S=?*J-1 

SUMSSaSUMSS+lDO/ ( S**?*(R**?+ ( p*S )**?)** ? ) 

1  s  SUMS  =  SUMS+1  DO  /  (  (  R**?  +  (  0*9  )  **?  )  **7  ) 

F®  T-l 
R  =  2*  T-l 
S= 2*1-1 

A ( I  ♦  K  )  =  1  DO-SUMP 
A ( K 9  I ) =SUMS*P**4 

A(!,M)=4PO0/(PT**4) *P**2* ( -1  DO ) **F / ( S**2 )* ( 1 D0-?4D0*SUMRR / ( PI**? )  ) 
A  ( t..  *  M  )  =  0  D  0 
N=  I  -1 

A(K»M)=-11 52D0/ ( P I **6 ) *P** 2* SUMSS* ( -1  DO ) **E 
2  CONTINUE 
N  =  ?*L 
F  =  L-1 

A ( L *N ) = ( -IDO ) **E 
M=?*L+1 
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T  T  =6 

CALL  JORDAN  (N*M*A*IND*TT  ) 

IF ( IND.FQ.O )GO  TO  66 
WR I TF ( 6  *  1 0 ) 

10  PORMAT <24H  THF  MATRIX  I S  SINGULAR.) 

STOP 

66  DO  60  T  =  1  * L 
K  =  L+I 

R  (  I  )  =  A  (  I  *  M ) 

6 D  C(n=A(K*M) 

6  WR I TF ( 6  *  1 1  )  P 

11  FORMAT ( 1H1 *  3H  P  =  Dfl.2*27H  B(R)  ) 

WRITF(6*1?)  (R(K  )  *K  =  1 *L> 

12  FORMAT (  1H  *  18X*D24.16) 

E  =  L  +  1 

WRITE (6 *33) 

33  FORMAT (1HK*3RH  C(S)  ) 

WRITE (6*1?) ( C ( K ) *X=F  *  N ) 

DO  16  1  =  1  *  L 
R=2* 1-1 
F=  I  —  1 

D017  J= 1 *L 
S=2* J-l 
X  =  J+L 
F  =  J-l 

17  ARS  (  I  *  J  )  =(-P**4*  (-1  DO  )  **F*R  (  I)  +  (  -1  DO)  **F*(  (P**?  +  -anO*  (  P*S  )  **?  )  *C  (K  ) 
1  -P**2*llS?D0/(Pt**6 )* (-1D0 ) **F/ ( S**? ) ) )/( ( R**2  + ( P*S ) **2 ) **2 ) 

16  CONTINUE 

36  FORMAT ( 1HK  *5  8H  SUMREQ 

1SUMSE0) 

WRITE (6 *36) 

DO  74  I  =  1  ♦  L 
c:mmr  =  oDo 
C(JMS  =  0D0 
F=  T-l 
K=T+L 
6  =  2  *  I  - 1 
DO  36  J=1 *L 
F  =  J-l 
R=2* J-l 

SUMR  =SUMR+(-lD0)**F*(R**2-(P*S)**2 )*ARS( J*T ) 

7R  6UMS  =  SUMS+(-lD0 )**E*ARS ( I  * J ) 

SUMR=-SUMR+C( K )  -4RD0/ (PT**4 ) *P**2* (-1 DO)**F/ ( S**2 ) 

37  FORMAT ( 1H,  1RX*2D24.16) 

74  WRITE (6 *37)  SUMR* SUMS 

FUMC=OD0 

no  3  P  1  =  1  *L 
F*  T-l 
X  =  L+  I 

79  FUMS=SUMS+C ( K ) * ( -1  DO ) **F 
WRITE (6 *40) 

40  FOPMAT  (1  HK  *23H  SIJMC  (  S  )  **  (  S-l  )  /?  ) 

WRITE (6 *37)  SUMS 
WR  I T  F ( 6  *  1 3 )  P 

13  FORMAT ( 1 H3  *  3H  P=  F6.2*^RH  THF  COFFF I C  TENTS  ARS  C(S).NE.O) 
00  14  1=1*5 

26  FORMAT ( 1HK  *6F1? .6 ) 
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14  WRITE (6 *26) (  ARS( T  * J>  * J=1  ,3  > 

C  CALCULATION  OF  THF  X  -  VELOCITY  COMPONENT 

WRITE (6 *22)  P 

22  EORMAT  ( 1 H2  * 3H  P=  Ef,.?,??H  X  -  VELOCITY  ) 

WR ITF ( 6  *  28 ) 

2 8  EORmaT(1HX*  7 OH  Y  =  0  PI/20  2PT/20  7PT/70  4PI/20  SPJ/20  6PI/20  7 

1PI/20  8PI/20  9PI/20) 

X  =  — P I  /I ODO 
D0 1 R  N  =  1 *6 
E*N-1 

WRITE (6 *29)  E 

70  FORMAT (1HK*2HX=I1*5HPI /10) 

X=X+P I / 1 ODO 
Y=— P I /?0D0 
0019  X  =  1  *10 
Y=Y+P! /20D0 
U$(X )=OD0 
DO 20  1=1 *L 
UR (K ) =000 
D021  J=1 *L 
R=2*J-1 

I F  ( N.LT  #6 ) GO  TO  74 
E  =  J— 1 

UR  (X  )  =UR(!C)+ARS(  J*  I  )*  (-100>**f 
GO  TO  21 

24  UR ( K ) =UR ( K ) +ARS ( J  *  I >*SIN(R*X) 

21  CONTINUE 
S=2*I-1 

70  US (X )=US (K)+UR (K)*S*rOS(S*Y ) 

UP ( X )=600*p*x*( (Y/PI ) **2-.2SD0) /PI 
19  U  (  X  )  =  PI/2D0*P*IJS(K) 

WRITF(6*2S) (U(X  )  *  X  =  1  *10) 

27  FORMAT! 1HJ*10D13. 5 ) 

28  FORMAT! 1HJ,10F7.^) 

WRITE (6 *30) 

30  FORMAT ( 1 HJ *40H  PARABOLIC  ) 

IP  WRITE(6*25) <UP(X) ,X*1 *10> 

C  CALCULATION  OF  THE  Y  -  VELOCITY  COMPONENT 

WR  TTE ( R  *  3 1  ) 

31  FORMAT ( 1H2 *SSH  Y  -  VELOCITY  > 

WR ITE ( 6 ♦ R2 ) 

82  EORMAT (  1HX  *  70H  X  =  0  PI/20  2PI/20  3PI/20  4PI/20  5PI/20  6PI/20  7 

1PI/20  8PI/20  9PI/20) 

Y=— P I / 1 ODO 
DO  78  N= 1 *6 
E  =  N-1 

WRITE (6 *79)  E 

79  EORMAT ( 1HK  *  2HY  =11 *  5HP  I / 1 0 ) 

Y  =  Y+P  T  /I  ODO 
X=-P I /20D0 
DO  80  X  =  1 *10 
X=X+P I /20D0 
VR (X ) =OD0 
DO  70  1  =  1, L 
VS (X ) =000 
no  71  J  =  1  ,L 
S=7* J—l 
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71  VS(K  )=VS(K)+ARS(  T  *j)*S!N(S*Y) 

P=7*T-1 

70  VR ( x ) =VP ( K ) +VS ( K ) *R*C0S ( R*X > 

VP (X )=S00*Y*( ,?^nn-( Y/PI )**?/70O)/PT 
P0  V  ( K  ) =-P I /?DO*VR ( K ) 

WR I TF ( 6  *  25 ) (V(K) *K  =  1 *10) 

WR I TF ( 6  * 30 ) 

7 P  WRITF(6*?5)  (VP( K)  ,K  =  1  *10) 

Y=-d I /I 0D0 
00  P5  N  =  1  *6 
v=v+dt /i ono 
V(N)=ODO 
no  50  T  =  1  *L 
K  =  L+ 1 
S = ?  *  I  - 1 

Ro  V(  M  )  =  V (  N  )  +C  (  K  )  IN'  (  v  ) 

P 5  FONT T MUR 
WRTTF(5*51 ) 

5 1  FORMAT ( 1 HK  *  24H  Y-VFLOCITY  AT  X=PT/2) 
WR ITF ( 6 ♦ ?P ) ( V ( M ) *M=1 *5) 

100  CONTIMtlF 
7  STOP 

fnd 
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j  #  r  •=■  t  o*  on 
T  A  IsM 
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(  V  -"l/T  -  »'  (  ,p4|l^w.(mw  O  a 
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f  .  M  ■’TT  «w 

* v  t  ■  -  y  r  fl  v/TT'-;  .-V  '  *H  f ) T AMflOl  r  a 
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SUBROUTINE  JORDAN ( N *M , A  *  I  NO , IT) 


C 

C 

C 

C 

c 

c 


DOUBLE  PRECISION  A(50»5lNP»B 
INTFGER  LR ( 50 ) *LC(51 ) 

M  IS  NO.  OF  COLUMNS 
N  IS  NO.  OF  ROWS 
A  IS  MATRIX 

IF  IND  IS  ZERO*  INVERSION  IS  GIVFN 
IF  IND  IS  ONE*  MATRIX  IS  SINGULAR 
IT  TS  OUTPUT  UNIT 


3  KZ  =  1 


DO  16  K= 1 *N 
LR (K )*K 
LC ( K ) =K 
T  =K 
I Z  =  K 

P  =  DABS( A (K*K)  ) 

DO  5  J  =  K*N 
DO  5  JZ  =K  ♦  N 

I  F ( P .GF • DABS ( A ( J  *  JZ )  )  ) GO  TO  5 
P  =  DABS ( A ( J  *  JZ )  ) 

I=J 

IZ=JZ 

5  CONTINUE 

I F ( I .GT • K ) GO  TO  23 

6  IF ( I Z.GT .K ) GO  TO  IB 

7  lF(P.GF.l.onoooonnonnnnnnD-?4)GO  to  i? 
WRITF(IT*42) 

I  ND  =  1 
RF  TURN 

12  P  =  i  .onnooonnnoooori^D^ /A  ( K  *K  ) 
A(K^)sl.on^nnnnnnnnnonoDn 
DO  13  J=1*M 

13  A(K*J)=P*A(K*J) 

DO  16  I =1 *N 

IF (X .FO. I ) GO  TO  16 

14  R= A ( I *K ) 

A  (  I  ,x)=n.noononnooonr»onODo 
DO  15  J= 1 *M 

15  A( I *J)=A ( I *J)-R*A(K*J) 

16  CONTINUE 
KZ  =  0 

X  =  N 

DO  26  I  A  =  1 *N 

IF (LR(K  > .EQ.K ) GO  TO  21 

17  I Z  =LR ( K  ) 

C  CHANGE  COLUMNS  IZ  AND  K 

18  DO  19  J= 1 *N 
B=A( J*  TZ  ) 

A ( J* IZ ) =A ( J*K  ) 

19  A  (  J  *  K  ) =B 
IF(KZ.FO.O) GO  TO  21 

20  LC ( K ) = I Z 
GO  TO  12 

21  IF (LC(K ) .FQ.K ) GO  TO  26 

22  I =LC ( K ) 
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CHANGE  ROWS  I 

AND  K 

23 

DO  ?4  J  =  1 *M 

R=A  (  I  * J ) 

A (  I  * J ) = A (K,J) 

24 

A( K * J)=R 

IF(KZ.FQ.O)GO 

TO  26 

25 

LR(K ) =1 

GO  TO  6 

26 

K=K-1 

IND  =  0 

RETURN 

42 

FORMAT ( 24H  THE 
END 

MATRIX  IS  SINGULAR 

i 


rc 


ot 


(  .C  •  ••  l.-U*  ’  ?  ^  T  *!'-TT  AM 


t  -,ci 

M.  r  =! 

(I  ♦  :  )  A  S  a 
(  '  •  >)  A  =  f '  .  T  >  ' 

v  )  a  4i  c 

.  .  '■  N  )  ~  T 

T  =  <  > )  «  J  p  c 

*  o  t  on 

r-N=N  r5  s 

o=ni4T 
‘•'C'  T^q 
-jut  h^CITAMOO^ 

01*? 
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P=  1.00 

THE  COEFFICIENTS  AKS 

CIS) .NE.O 

-0.016156 

-0.011838 

0.003941 

0. 000303 

0.000004 

0.006189 

0.003445 

-0.002  }84 

-0.000273 

0.000008 

-0.002461 

-0.C00904 

0.001201 

0. 000215 

-0.000016 

• 

0.001271 

0.000294 

-0.000626 

-0. 000158 

t 

0.000018 

-0.000779 

-0. 0001 08 

0.000351 

0.000119 

-0  .000020 

0 


, 


i 


*  ■ 


■' . ' 


!  ■ 


f  * 


! 


• 
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P= _ 1.00 _ X  -  VELD C ITY 


Y=0  PI/20  2PI/20  3PI/20  4PI/20  5PI/20  6PI/20  7PI/20  8PI/20  9PI/20 


X=0PI/10 

t  , 

0.000  -0.000  0.000  -0.000 

0.000 

- •  tfl - 

-0.000 

0.000  -0.000 

0.000 

-0.000 

PARABOLIC 

-0.000  -0.000  -0.000  -0.000 

-0.000 

-0.000 

-0.000  -o.coo 

-0.000 

-0.000 

X=1PI/10 


-0.009 

-0.009  -0.007  -0.004 

-0.001 

0.003 

0.006 

0.C07 

0.006 

0.004 

PARABOLIC 

-0.150 

-0.148  -0.144  -0.136 

-0.126 

-0.112 

-0.096 

-0.076 

-C. 054 

-0.028 

X  =  2P  I  /  10 


-0.022  -0.021  -0.017  -0.011 

-0.003 

0.006 

0.013 

0.017 

0.016 

0.010 

PARABOLIC 

-0.300  -0.297  -0.288  -0.273 

-0.252 

-0.225 

-0.192 

-0. 153 

-0.108 

-0.057 

X  =  3P  I  / 10 

- .  .  /  -- 

-0.039  -0.039  -C.034  -0.027 

> 

-0.011 

0.007 

0.026 

0.C36 

0.036 

0.022 

PARABOL  IC 

-0.450  -0.445  -0.432  -0.409 

-0.378 

-0.337 

-0.288 

-0.229 

-0.162 

-0.085 

I 

X  =  4P  I  /  1 0 

-0.052  -0.056  -0.063  -0.059 

-0.037 

0.002 

0.044 

0.071 

0.073 

0.046 

PARABOLIC 

-0.600  -0.594  -0.576  -0.546 

-0.504 

-0.450 

-0.384 

-0.306 

-0.216 

-0.114 

*  , 

X  =  5P  I  /  10 

i 

■- 

-0.167  0.070  -0.219  -0.001 

-0.194 

0.074 

-0.012 

0.171 

0.068 

0.083 

PARABOLIC 

-0.750  -0.742  -0.720  -0.682 

-0.630 

-0.562 

-0.480 

-0.382 

-0.270 

-0. 142 

' 


. 


•  •  o  ■■  ’  n:.m 


o.  = 


■ 


.  -  vis.o-  oso.<mxJ^m 
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Y  -  VELOCITY 


x=o 

PI/20 

2PI/20 

3P I / 20 

4P  1/20 

5P  1/20 

6PI /20 

7P I /2  0 

8  P  I  /  2  0 

9P  I / 20 

Y  =  OP  I  /  10 

-0.000 

0.000 

-o.coo 

0.000 

-0.000 

0.000 

-0.000 

0.000 

-0.000 

0.000 

PARABOL  IC 

0.000 

0.000 

o.coo 

0.000 

0.000 

0.000 

0.000 

o.coo 

0.000 

0.000 

Y=1PI/10 


0.009 

0.007  0.011  0.010 

0.016 

0.016 

0.022 

0.017 

0.019 

-0.011 

PARABOLIC 

0.148 

0.148  0.148  0.148 

0 . 1 48 

0.148 

0.14  8 

0.  148 

0.  148 

0.  148 

Y=2PI/10 


0.011 

0.012  0.014  0.018 

0.024 

0.032 

0.041 

0.050 

0.058 

0.051 

PARABOL  IC 

0.284 

0.284  0.284  0.284 

0 .284 

0.284 

0.284 

0.284 

0.284 

0.284 

Y=3PI/10 


0.009 

0.010  0.011  0.016 

0.021 

0.030 

0.040 

0.055 

0.069 

0.070 

PARABOLIC 

0.396 

0.396  0.396  0.396 

0.396 

0.396 

0.396 

0.  396 

0.396 

0.396 

Y  =  4  P I /  10 


0.004 

0.002  0.005  0.004 

0.009 

0.010 

0.018 

0.019 

0.033 

0.018 

0.472 

PARABOL  IC 

0.472  0.472  0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

Y  =  5P  I  / 10 

- 

0.000 

0.000  O.COO  0.000 

O.COO 

0.000 

0.000 

0.000 

0.000 

0.000 

0.500 

PARABOLIC 

0.500_ 0.500_ 0.500 

0.500 

0-500 

0-500 

Q-500 

Q-^on 

n^nn 

. 


> 
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FORTRAN  SOURCE  STATEMENT  FOR  THE  VELOCITY  FIELD  -  CS=0 

DOURLE  PRFCISION  A(50,51>*  P  I  ,P  ,R  ,S ,  SUMR  *  SUMRR  ,  SUMS  ,SIJMSS  *  PT  ,PTP , 
1  PR  (  10)  ♦  PRP  no  )  *  X » Y  *  W  »  Z  *  UR  (  1 0 ) ,US ( 1 0 ) *U < 1 0 ) * VR ( 10 ) , VS ( 1 0  > , V < 1 0  )  * 

1  UPU0)*VP(l0),C(?5>,RI(?5)  *BI I ( 25  >  *AI ( 25  >  *ARS(50*50  >  *<  I ( ?5 ) ,R(25  > 
INTEGER  I ,J,K,L ,M,N ,RHO,F,F ,H 
PI=3. 141592653589793 
DO  100  RHO= 1 ,5,2 
P  =  RHO 
L  =  1 0 

DO?  I  =  1  *  L 
K  =  L+  I 
D03  J=1 *  L 
E=  I + J-2 
R=?*J-1 
5=?*I-l 

A ( I  * J)  =  ( -IDO )**E* ( R**?-(P*S )**?)/(  ( R**2  + ( P*S )**?)**? )*P**4 
R=?*I-1 
S=?*J-1 
N  =  L+J 

A ( K ,N ) =- (-] DO ) **F* ( R* *2+3 DO* ( P*S )**?)/(  ( R**?  + (P*S )**?)**?) 

A ( I ,N ) =0D0 
3  A  (  K  ,  J  ) =  ODO 
SUMRsODO 
SUMRR=ODO 
$UMS=ODO 
SUMSS=ODO 
DO  15  J=1 ,L 
S=?*I-1 
R=?*J-1 

SUMRR=SUMRR+(R**2-(P*S)**? ) /( <R**?+(P*S>**2>**2> 

SUMR=SUMR+( (R**?— (P*S)**?)*(R**2+3D0*(P*S)**2 ) ) /( ( R**2+ ( P*S ) **? ) ** 
1?  ) 

R  =  2* I  —  1 
S  =  2* J-l 

SUMSS=SUMSS+1 DO/  ( S**?* ( R**?+ ( P*S )**?)**? ) 

15  SUMS=SUMS+1 DO / ( ( R**?  + ( P*S > **2 ) **2 ) 

F  =  I  —  1 
R=2* I -1 
S  =  2  *  I — 1 

A ( I  *  K ) - 1  DO— SUMR 
A  (  K  *  I ) =SUMS*P**4 
M=?*L+1 

A (  T ,M)=48D0/( PI**4 )*P**?*(-1 DO ) **F / ( S**2 ) * ( 1 DO-?4DO*^UI  RR/(PI**2) ) 
N=  1-1 

A ( K*M)=-115?D0/ ( PI**6 )*P**2*SUMSS*(-1D0 ) **E 
2  CONTINUE 
N  =  ?*L 

DO  50  1=1, L 
K  =  L+ 1 

R ( I )=A(K,M)/A(K,I  ) 

50  C  ( K  ) =0D0 
WR I TE ( 6 , 53 ) 

53  FORMAT ( 1HK,7?H  SUMR  48D0/IPI 

1**4*(-1D0)**E/(S**?>  > 

DO  51  I =1 *L 


r  '*JH  v  -  T  |  )v  ih  T 


.  3  M  T  •  r  r.  :  ■>  <r,»  I  '•  '  1/  :  9  T  qo  9 


,nro,  ro,  a  ,  *>  .  .  .  ,  ’  .  {  r  =>  ,  r  ~  )  A  ‘  !  T  ' '  r  ■  1  °  n  I  f  r !  1  o  n 

#  /  r  )  -  •  .  vr-'/»  n  r  -  •  '  .  (  o  r  .  *  •  f  •  ->  .  (  r  >  .  .  .  Y  .  >’  *  (  r  )  :.qa,  0  r  )  q<  f 

.  -  )  '  T  A  .  -  .  '  r  |  .  .  f )  qtJ  I 

.  .  i.  u  r  .  #  .  .  ,  .  '  ;:v.  r  m 

r  \  \  .  1  p  ,r- 

,<>,  fsouc  r  r  on 
0HO  =  q 
n  r=  ! 
jt r=T  ?oo 
T  +  ]  =  * 
j. r=t  *oo 
9- L«M  =1 

,>*  fq*(r*#(rM(?#q)+r  #q  )  )  \  (  Ct**  (  3*q  )  -C,  *  *  \  ^  i  *  •  { on  T ~ )  =  f  L *  T  )  A 

r_T#<:  =  cj 

t+ 

,  -w.  '  r~  v *  (  -. * n  )  \  (c#*f  =  a 

0 .1  o  -  (  ^  #  t  )  r 
ono=(t#^)A  r 
ono=qvi \-> 
nop  =  qflMt  *5 
Of!0  =  »Mf !? 

J*  f *L  3  f  on 

r-L*<;=n 

'  -  *  ■■-  f  c  ■■  v-  {  -  ■'  a  >  4-  *  '^n\(  r'5f  (  3  ••*  1  )  -c  •  *  '  )  +<3--:  'M  |-»-qqM»p 
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r-i 
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r-!=q 
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K  =  L+  I 
F=  1-1 
S=2*I-1 
St  IMR-ODO 
DO  52  J  =  1  *L 
R=2* J-l 

5?  SUMR=SUMR+A ( I  * J ) *R ( J )  +  ( -IDO ) **F*1152DO/ ( PI**6 ) *( R**?-( P*s )**?)/( (R 
1**2  + (p*s )**? )**2 ) / ( S**2 ) 

SUMS=48D0/(P 1**4 )*P**2*( -IDO )**£/( 5**? ) 

51  WRITE (6 *37)  SUMR  *  FUMS 
6  WRITE (6 *11)  P 

11  FORMAT ( 1H1 *3H  P=  D8#2*27H  BCR)  ) 

WRITF(6*12) ( R ( K ) *K  =  1 *L> 

12  FORMAT ( 1H *  18X*D24.16) 

F  =  L  +  1 

WRITE (6*33) 

33  FORMAT ( 1HK  *39H  C(S>  ) 

WR I TF ( 6  *  1 2 )  (C(K)  *K=F*N> 

DO  1. 6  1  =  1  *L 
R=2* 1-1 
F=T-1 

DO 1 7  J  =  1  ♦  L 
5=2* J-l 
K  =  J+L 
F  =  J-l 

17  A  R  S (  I  *  J )  =  (-P**4* (-1  DO ) **F*R (  I)  +  ( -1  DO ) **F* (  ( R**p+3D0* ( P*S ) **2 ) *C( K ) 
1  — P**2*1152D0/(PI**6)*(-1D0)**F/(S**2) ) )/( ( R**2  + ( P*S ) **2 ) **2 ) 

16  CONTINUE 

36  FORMAT ( 1 HK  *  5  8H  SUMRFQ 

1SUMSEO) 

WR  ITF ( 6  * 36 ) 

DO  34  I  =  1  *L 
5UMR=0D0 
5UMS=0D0 
F=  T-l 
K=I+L 
S=2*I-1 
DO  S5  J= 1  * L 
F  =  J-l 
R  =  2  *  J— 1 

SUMR  =SUMR+ ( -IDO ) **F* ( R**2" ( P*S ) **? ) * ARS ( J , I ) 

35  SUMS = SUM S+ ( -1  DO ) **F*AR5 ( I  * J) 

SUMR=— SUMR+C ( K )  -48DO/ (P 1**4 ) *P**2* (-1D0 >**F/ ( S**2 ) 

37  FORMAT ( 1H*  18X*2D24.16> 

?4  WR I TF ( 6  *  37 )  SUMR* SUMS 

WR ITF ( 6  *  13 )  P 

13  FORMAT ( 1H3  *  3H  P=  F6.2*38H  THE  COEFFICIENTS  ARS  CCS).  F.O) 
DO  14  1=1*5 

26  FORMAT ( 1HK*5F12*6) 

14  WRITF(6*26) ( ARS( I ♦ J) * J=1 *5 ) 

C  CALCULATION  OF  THE  X  -  VELOCITY  COMPONENT 

WR  I TE ( 6  *  22 )  P 

22  FORMAT ( 1 H2  *  3H  P=  F6.2*23H  X  -  VELOCITY  ) 

WR  ITF ( 6  *  28 ) 

28  FORMAT ( 1 HK  *  70H  Y  =  0  PI/20  2PI/20  3PI/20  4PI/20  5PI/20  6^1/20  7 

1PI/20  8PI/20  9PI/20) 

X=-PI/10D0 


U  )=N 

r_T  -  "> 
o  q  °  -  o  v  •  i  “ 

u  =t  n 

f-|.*C  =  a 

- ;  )  )  V  (  ■  ■  *  t  ->  J.  q  )  _  •  i  •?  .  \  I'l'-i  r  r^^!'  (  r  ’  -  '  4-  /  ;  \  ■  r’  (  !  >  T  •  A  *■ .  1  '  ?  -  ’  !  °  ' ? 

(Cy*'p  I  *  *  (  c  #  M  p  *P  )+<*,**  T 

(  y  •<  t  )  \  ,  v  •  (  a  r  r  -  )  *  -  *  x  i  '  '  ■  -r  O  )  \  >  n  p  -  '  •  1  ? 

#  n (  rc  .  _>  >  “IT  j  PM  r, a 

n  /  r  .  >  *  TT  P'-'  A 

(  i-iq  JV  ,  c#pf|  -a  H p  •  I H  r  )  T  A  M9101  rr 

(  i .  r  -  * .  f  n  >  <  )  (  c  .  *  >  -3T  t  Pl ' 

(  s  f  ,Ar  i * x h  r  .  -« r  )  rAMfln^  c  r 

( ' . *  )1T|  qw 

{  hpp  »  xh  r )  t  AMqm  pp 

(  »p  !=  N*  (  •  )  '  (  C  !  .  '  "  T  T  c  .• 

j •  r  =  T  \  r  or 

f  —  T *Crq 

r-7  =q 

j#  r*t  v 
I~L*S*£ 
HL** 
r-t*q 
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,  r  .  I  )  ->  Q  a  .  r  *•».'.">  *  Q  )  -  c  Q  )  .  v-  (  or\  r  ~  )  +<-  ip=  vMi)? 
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0018  N=l*6 
E  =  N-1 

WR I TE ! 6 ♦ 29 )  E 

29  FORMAT !1HK*2HX=I1,5HPI/10) 

X  =  X  +  P I / 1 ODO 

Y  =  *-P  I  /20D0 
DO 1 9  K=l *10 

Y  =  Y+P 1/2  ODO 
t  IE  (K  )  sODO 
D020  T  =  1  *L 
UR ( K ) =0D0 
DO? 1  J=1  *  L 
R=2* J-l 

21  UR ( K ) =UR ( K ) +ARS ( J  ,  I >*SIN(R*X) 

S  =  2  *  I  - 1 

20  US ( K ) =US ( K ) +UR (lO*S*C0S(S*Y ) 

UP  ( K  ) =6D0*P*X* ( ( Y/PI ) **2  — • 2  EDO ) /PI 
19  U ( K ) =  PI/2D0*P*US(K> 

WR  I TE  (6*25)  !U(K)  *K  =  1  *10) 

23  FORMAT! 1HJ,10D13. 5 > 

?5  FORMAT! 1HJ*10F7. 8) 

WRITE!*  *30) 

30  FORMAT!  1HJ,/+0H  PARABOLIC  ) 

18  WRITE(6*25> !UP(K> *K=1*10> 

C  CALCULATION  OF  THE  Y  -  VELOCITY  COMPONENT 

WR I TE ( 6  *  31 ) 

31  FORMAT! 1H2*35H  Y  -  VELOCITY  > 

WR I TE ( 6  *  82  ) 

82  FORMAT (  1 HK  *  70H  X  =  0  PI/20  2PI/20  3PI/20  4PI/20  5PI/20  6PI/20  7 

1PI/20  8PI/20  9PI/20) 

Y  =  — P I /I ODO 
DO  78  N  =  1  * 6 
E  =  N-1 

WR I TE ( *  *  79 )  E 

79  FORMAT ! 1HK*2HY=I 1 *  5HP  I / 1 0  ) 

Y=Y+P I / 1 ODO 

X  =— P I /20D0 
DO  80  K  =  1 ,10 
X  =  X  +P  1/2  ODO 
VR!K)=ODO 
DO  70  I  =  1  *L 
VS ( K ) «0D0 
DO  71  J  =  1  ,L 
S=2* J-l 

71  VS ( K ) =VS ( K ) +ARS ( I *J>*SIN(S*Y) 

R  =  2* I -1 

70  VR(K>=VR!K>+VS(K>*R*COS(R*X > 

VP(K)=6D0*Y*(.25D0-!Y/PI )**2/3Dn) /PI 

80  V ( K ) =— P I /2D0*VR ! K ) 

WRITE! 6 *25) (V(K> *K=1*10> 

WR I TR ! 6 , 30 ) 

78  WRITE(6*25>  ! VP ( K  > *K  =  1 *10> 

100  CONTINUE 
7  STOP 
END 
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P=  1.00  THE  COEFFICIENTS  ARS  C(S).  E.O 


-0.012555  -0.010149  0.001627  -0.000449  0.000169 


0.003400  0.002238  -0.000701  0.000248  -0.000104 


-0.000866  -0.C00415  0.000226  -0.000107  0.000053 


0.000291  0.C00101  -0.000077  0.000047  -0.000027 


-0.000121  -0.0C0031  0.0C0030  -0.000021  0.000014 


' 


. 


V 
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P=  1.00  X  -  VELOCITY 


Y  =  0  PI/20  2PI/20  3PI/20  4PI/20 

5P  1/20 

6PI/20  7PI/20 

8PI/20  9PI/20 

X  =  OP  I / 10 

0.000  -0.000  0.000  -0.000  0.000 

-0.000 

0.000  -0.000 

0.000  -0.000 

PARABOLIC 

-0.000  -0.000  -O.CCO  -0.000  -0.000 

-0.000 

-0.000  -0.000 

-0.000  -0.000 

X=1PI/10 

-0.010  -0.009  -0.007  -0.004  -0.001 

0.003 

0.006  0.007 

0.006  0.004 

PARABOL IC 

-0.150  -0.148  -0.144  -0.136  -0.126 

-0.112 

-0.096  -0.076 

-0.054  -0.028 

X=2P I / 10 

C 

-0.023  -0.021  -0.017  -0.010  -0.002 

0.006 

0.013  0. 017 

0.016  0.010 

PARABOLIC 

-0.300  -0.297  -0.288  -0.273  -0.252 

-0.225 

-0.192  -0.153 

-0.108  -0.057 

X=3P I / 10 

-0.041  -0.039  -0.031  -0.020  -0.006 

0.009 

0.022  0.031 

0.031  0.021 

PARABOLIC 

-0.450  -0.445  -0.432  -0.400  -0.378  -0.337  -0.288  -0.229  -0.162  -0.085 


X  =4P I / 10 


-0.062  -0.058  -0.049  -0.033 

-0.013 

0.009 

0.030 

0.047 

0.053 

0.040 

PARABOLIC 

-0.600  -0.594  -0.576  -0.546 

-0.504 

-0.450 

-0.384 

-0. 306 

-0.216 

-0.114 

X=5P I / 10 

— 

-0.074  -0.070  -0.059  -0.041 

-0.020 

0.007 

0.032 

0.056 

0.066 

0.060 

PARABOLIC 

-0.750  -0.742  -0.720  -0.682 

-0.630 

-0 .562 

-0.480 

-0.382 

-0.270 

-0.142 

s 


'*■ 


r- 


•  . .  1 


0 


' 


= 
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Y  -  VELOCITY 


x=o 

PI/20  2PI/20  3PI/20 

4P I / 2  0 

5P  1/20 

6PI/20 

7PI  /20 

8PI/20 

9PI/20 

Y=0PI/10 

-0,000 

0.000  -0.000  0.000 

-0.000 

0.000 

-0.000 

0.000 

-0.000 

0.000 

0.000 

PARABOLIC 

0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Y= 1 P  I  /  10 

L 

0.C09 

0.009  0.010  0.012 

0.014 

0.017 

0.019 

0.020 

0.018 

0.013 

0. 148 

PARABOL  IC 

0.148  0.148  0.148 

0.148 

0.148 

0.148 

0.  148 

0.148 

0.148 

Y  =  2P  I  / 10 

0.012 

0.013  0.015  0.018 

0.022 

0.026 

0.030 

0.033 

0.031 

0.022 

0.284 

PARABOLIC 

0.284  0.284  0.284 

0.284 

0.284 

0.284 

0.284 

0.284 

0.284 

Y  =  3P  I  /  10 

0.010 

0.010  0.012  0.015 

> 

0.019 

0.024 

0.028 

0.033 

0.033 

0.025 

0.396 

PARABOLIC 

0.396  0.396  0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

Y=4PI/10 

0.003 

0.004  0.004  0.006 

0.007 

0.010 

0.013 

0.016 

0.019 

0.017 

0.472 

PARABOLIC 

0.472  0.472  0.472 

0.4  72 

0.472 

0.472 

0.472 

0.472 

0.472 

\ 

y  =  5p  I  /  10 

-0.000 

-0.000  -0.000  -0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

0.500 

PARABOL iC 

0.500  0.500  0.500 

0.500 

0.500 

0.500 

0.500 

0.500 

0.500 

■ 


P=  3.00 

THE  COEFFICIENTS  AUS 

C ( S) .  C.O 

-0.001691 

-0.001401 

0.000199 

-0.000053 

0.000020 

0.001395 

0.001128 

-0.000181 

0. 000050 

-0.000019 

-0.0C0977 

-0.000753 

0.000148 

-0. 000043 

0.000017 

0.000618 

0.C00444 

-0.0C01 10 

0. 000035 

-0.000014 

-0.000378 


-0 . C00249 


0.0C0078 


-0. 000028 


0.000012 


- 


.  > 
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P  =  3. CO  X  -  VELU C ITY 


Y  =  0 _ Pl/20  2P  1  /20  3P1/20  4PI/20  5PI/20  6PI/2Q  7PI/2Q  8PI/20  9P1/20 


X  =  OP  I  / 10 


0.000  -0.000  0.000  -0.000 

0.000 

-0.000 

0.000 

-0.000 

0.000 

-0.000 

PARABOLIC 

-0.000  -0.000  -0.000  -0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

X=1PI/10 

.M  .«•«*» 

0.000  0.000  0.000  0.000 

0.000 

0.000 

-0.000 

-0.000 

-0.000 

-0.000 

PARABUL IC 

-0.450  -0.4**5  -0.432  -0.409 

-0.378 

-0.337 

-0.288 

-0.229 

-0.162 

-0.085 

X  =  2P  I  / 10 

1 

c 

-o.coo  -0.000  -0.000  0.000 

0.000 

0.000 

0.000 

0.000 

-0.000 

-0.000 

PARABUL 1C 

-0.900  -0.B91  — 0 .864  -0.819 

-0.756 

-0.675 

-0.576 

-0.459 

-0.324 

-0.171 

X  =  3P  I  /  10 

-0.003  -0.002  -0.002  -0^001 

0.000 

0.001 

0.002 

0.002 

0.001 

0.001 

PARABUL IC 

-1.350  -1.336  -1.296  -1.228 

-1.134 

-1.012 

/ 

-0.864 

-0.688 

-0.486 

-0.256 

X=4PI/10 

-0.024  -0.022  -0.018  -0.011 

-0.002 

0.007 

0.014 

0.017 

0.016 

0.010 

PARABOL IC 

-1.800  -1.782  -1.728  -1.638 

-1.512 

-1 . 350 

-1.152 

-0.918 

-0.648 

-0.342 

X=5PI/10 

\ 

- 

-0.073  -0.069  -0.059  -0.041 

-0.019 

0.008 

0.033 

0.056 

0.066 

0.056 

PARABUL IC 

-2.250  -2.227  -2.160  -2.047 

-1.890 

-1.687 

-1.440 

-1.147 

-0.810 

-0.427 
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Y  -  VELOCITY 


x=o 

PI/20  2PI/20  3PI/20 

4PI/20 

5PI/20 

6PI/20 

7P I /2  0 

8PI/20 

9PI/20 

Y  =  0  P  I  /  10 

-0.000 

0.000  -0.C00  0.000 

-0.000 

0.000 

-0.000 

0.000 

-0.000 

0.000 

0.000 

PARABOL IC 

0.000  0.CG0  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Y  =  1 P  I  /  10 

-0.000 

0.000  -0.000  0.000 

-0.000 

0.001 

0.002 

0.006 

0.013 

0.021 

0.148 

PARABOL  IC 

0.148  0.148  0.148 

0. 148 

0.148 

0.148 

0.  148 

0.148 

0.148 

Y  =  2P  I  /  10 

-0.001 

0.000  -0.001  0.001 

-0.000 

0.002 

0.002 

0.010 

0.019 

0.035 

0.284 

PARABOLIC 

0.284  0.284  0.284 

0.2  84 

0.284 

0.284 

0.284 

0.284 

0.284 

Y=3P  I  /  10 

• 

i 

pH 

o 

o 

• 

o 

1 

0.001  -0.001  0.001 

-0.001 

0.002 

0.001 

0.008 

0.016 

0.035 

0.396 

PARABOL  IC 

0.396  0.396  0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

Y  =  4P  I  /  10 

-0.001 

0.001  -0.001  0.001 

-0.001 

0.001 

-0.000 

0.004 

0.005 

0.019 

0.472 

PARABOLIC 

0.472  0.472  0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

Y  =  5P  I  / 10 

i 

— 

0.000 

0.000  0.C00  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.500 

parAbOl id 

0.500  0.500  0.500 

0.500 

0.500 

0.500 

O'.  500 

0.500 

0.500 
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P=  5.00 

THE  COEFFICIENTS  Arts 

C  (  S  )  .  E  .  0 

-0.000619 

-0.C00513 

0.0C0072 

-0. 000019 

0.000007 

0.000576 

0.C00474 

-0.0C0070 

0. 000019 

-0.000007 

-0.0C0502 

-0.C00406 

0.0C0065 

-0. 000018 

0.000007 

0.000412 

0.C00325 

-0.0C0058 

0.000017 

-0.000006 

-0.000323 


-0.000246 


0.0C0050 


-0. 000015 


0.000006 


1 


*  -  . 


57 


P=  5.00 

X  - 

VELOCITY 

Y  =  0 

PI/20 

2PI/20  3PI/20 

4PI/20  5PI/20  6PI/20  7PI/20 

8PI/20  9PI/20 

X=0PI/10 

0.000  - 

0.000 

O.CGO  -0.000 

0.000  -0.000  0.000  -0.000 

0.000  -0.000 

-0.000  - 

0.000 

PARABUL  IC 

-0.000  -0.000 

-0.000  -0.000  -0.000  -0.000 

-0.000  -0.000 

X=1PI/10 

0.000 

0.000 

O.COO  0.000 

0.000  0.000  -0.000  -0.000 

-0.000  -0.000 

-0.750  - 

0.742 

PARABOLIC 

-0.720  -0.682 

-0.630  -0.563  -0.480  -0.382 

-0.270  -0.142 

X=2PI/10 

s 

0.000 

0.000 

O.COO  0.000 

0.000  0.000  -0.000  -0.000 

l  , 

-0.000  -0.000 

PARABOLIC 

-1.500  -1.485  -1.440  -1.365  -1.260  -1.125  -0.960  -0.765  -0.540  -0.285 


0.000  0.000  0.000  0.000 

> 

0.000 

0.000 

-0.000 

-0. 000 

-0.000 

-0.000 

PARABOLIC 

-2.250  -2.227  -2.160  -2.047 

-1.890 

-1.688 

-1.440 

-1.147 

-0.810 

-0.427 

X  =4P I / 10 

-0.006  -0.005  -0.004  -0.002 

0.000 

0.002 

0.003 

0.004 

0.003 

0.001 

PAR  A BO L IC 

-3.000  -2.970  -2.880  -2.730 

-2.520 

-2.250 

-1.920 

-1.530 

o 

CO 

o 

• 

1 

-0.570 

X  =  5P  I  / 10 

l 

- 

-0.070  -0.066  -0.056  -0.038 

-0.017 

• 

0.009 

0.033 

0.054 

0.061 

0.050 

PARABOLIC 

-3.750  -3.712  -3.600  -3.412 

-3.150 

-2.812 

-2.400 

-1.912 

-1.350 

-0.712 

-  ;i  *  ~  r' '  *(lj 


■ 


.  4  • 


’ 


' 
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Y  ~  VELOCITY 


x=o 

PI/20  2PI/20  3PI/20 

4 P  I  / 20 

5P 1/20 

6P I / 2  0 

7  P I /2  0 

8  P I / 20 

9P1/20 

Y  =  OP  I  /  10 

-0.000 

0.000  -0.000  0.000 

-0.000 

0.000 

-0.000 

0.000 

-0.000 

0.000 

0.000 

PARABOLIC 

0.000  O.COO  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

( 

Y=1PI/10 

-0.001 

0.001  -0.001  0.001 

-0.001 

0.001 

-0.001 

0.002 

0.003 

0.019 

0.148 

PARABOLIC. 

0.148  0.148  0.140 

0.148 

0.148 

0.148 

0.  148 

0.  148 

0.148 

Y=2P I / 10 

-0.C01 

0.001  -0.001  0.001 

-0.001 

0.002 

-0.002 

0.003 

0.003 

0.030 

0.284 

PARABOLIC 

0.284  0.284  0.284 

0.284 

0.284 

0.284 

0.284 

0.284 

0.284 

Y  =  3P I /  10 

-0.001 

0.001  -0.001  0.001 

-0.002 

0.002 

-0.002 

0.003 

0.001 

0.029 

0.396 

PARABOLIC 

0.396  0.396  0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

0.396 

Y=4PI/10 

-0.001 

0.001  -0.001  0.001 

-0.001 

0.001 

-0.001 

0.002 

-0.001 

0.014 

0.472 

PARABOLIC 

0.472  0.472  0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

0.472 

Y=5PI/10 

- 

0.000 

0.000  O.COO  0.000 

0.000 

-0.000 

-0.000 

-0.000 

0.000 

0.000 

0.500 

PARABOLIC 

0.500  0.500  0.500 

0.500 

0.500 

0.500 

0.  500 

0.500 

0.500 

> 

. 


■ 


■ 
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